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Abstract 

We present the theory of an extremely correlated Fermi liquid with U — > oo. This liquid has an 
underlying auxiliary Fermi liquid Greens function that is further caparisoned by extreme correla- 
tions. The theory leads to two parallel hierarchies of equations that permit iterative approximations 
in a certain parameter. Preliminary results for the spectral functions display a broad background 
and a distinct T dependent left skew. An important energy scale A(fe, x) emerges as the average 
inelasticity of the FL Greens function, and influences the photoemission spectra profoundly. A 
duality is identified wherein a loss of coherence of the ECFL results from an excessively sharp FL. 



Introduction Correlated electron systems attract two distinct approaches. An interme- 
diate to strong coupling approach is used when the interaction U is comparable to the band 
width 2W, and has seen some success in recent timespQ. On the other hand, Anderson[2] has 
argued that myriad experiments on high T c superconductors require a better understand- 
ing of the t-J model physics. This model sets U — >■ oo right away i.e. leads to extreme 
correlations and involves Gutzwiller projected Fermi operators that are non canonical. 
Thus Wick's theorem is immediately lost, and perturbative schemes encoding the Feynman 
Dyson approach become useless. Since this approach is at the root of most current many 
body physics text books, the task of understanding the t-J model is not lightly undertaken. 

The Schwinger approach to interacting field theories is a powerful and attractive alterna- 
tive. It is fundamentally non perturbative, where Wick's theorem is bypassed by dealing with 
suitable inverse Greens functions. Conventional many body theory for canonical Fermions 
can also be cast into this approach, and leads to the standard results. In Ref. [3j (henceforth 
I), the author has recently applied the Schwinger method to the t-J model, and found a 
class of solutions that are termed as extremely correlated quantum liquids. That state is 
presumably realized under suitable conditions. However it gives a Fermi surface (FS) vol- 
ume that is always distinct from that of the Fermi gas. This is contrast to the case of Fermi 
liquids (FL), where the important theorem of Luttinger and Ward (L-W) [U [5] mandates 
the invariance of the FS volume under interactions. 

In this paper we propose a state of matter termed as an extremely correlated Fermi liquid 
(ECFL). The ECFL found here, represents an alternate class of solutions for the t-J model, 
where the Fermi surface satisfies the Fermi gas (i.e. L-W) volume. In this work we present 
the essentials of the formalism, and display preliminary results on spectral functions that are 
suggestive of the relevance of the ECFL state to cuprate materials. An inherent flexibility 
of the Schwinger approach permits the construction of an alternate class of solutions from 
the one found in I. The excitations of the ECFL state may be thought of as bare electrons 
undergoing a double layer of renormalization: the FL dressing into quasiparticles that are 
further caparisoned (i.e. decorated) by extreme correlations. 

Formalism: The physical projected electronic Greens function Q satisfies an equation 
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of motion (EOM) (1-29) written compactly in matrix form as 

(d Ti - n)g{i, f) = -8(i, f) {i - 7 «} - Vi ■ g(i, f) 
-x(i,]).g(if)-Y(ij).g(]j), (i) 

where /j, is the chemical potential and an implicit integration over space time variables such 
as j, written with bold overlined letters, is implied, 

X(i,j) = -t(i,j) (D{i) + D(j)) + ±J{i,Z) (D(i) + D(k))6(i,j) 

Y(i,j) = -t(i,j) (l- 7 (i)- 7 (j)) + Ij(i,k) (l- 7 (i)- 7 (k)){(i,j). (2) 

In the above expression^, we used 7 (i) = G k (i,i) with the k conjugation defined by 
(M k ) ai(72 = M^a x o\o%y and D ai(T2 {i) = <j\<Ji sv ^^ 2 ■ The added (Bosonic) source term 
Vf 1<J2 (Tj) is central to this approach; it is a space-time dependent field that couples to the 
charge and spin densities through a term in the action: ^ icr f Q dr Vf lfT2 (r) Xf icr2 (r), where 
Xp a2 is the spin and density operator at site i that acts as | cr^ ) (cr 2 1 . 

An important technical problem highlighted in I is to deal with the time dependence of 
the term in Eq. Q which makes the theory non canonical. Here we use the decompo- 
sition into two factors [7] : 

g(a,b) = g(a,b)-MM, (3) 

and express 7 (i) = (g(i,j) ■ fi(j,i)) k . The object g is an auxiliary FL Greens function and 
/i(b, b) is an appurtenant (or supplementary) factor that is determined below. Antiperiodic 
boundary conditions C?(0, 77) = —G(/3, rf) and Gfa, 0) = — G{t%, (3) imply that both factors g 
and n are Fourier transformed using Fermionic Matsubara frequencies. We define the inverse 
Greens function g _1 (a,b) • g(b, b) = 18(a,b), and thence a vertex function AJ^(p, q;r) = 
— sv $ 3 ct 4 {g^o-a (Pi o) }■ Thus g, fx and g _1 are matrices in the spin space, and the vertex A 
has four indices. We also define a linear operator 

Hi, f) = (t(i,J) r • g(j, /) - \j{hi) € ■ g(«, /)) + ( 4 ) 

where the matrix £* 1(T2 = <J\<Ji- The * is used as a place holder that transmits the spin indices 
(after conjugation) of the £ matrix to the source matrix V in the functional derivative. This 
notation used is illustrated in component form by ■ ■ ■ C a cr b ' ' ' = ■ ■ ■ a a ab ■ ■ ■ 5/5V? a " ah . 
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A useful chain rule for the functional derivative is noted 



D(r)-g(a,b) = f-g(a,c). A*(c,d;r) .0(d,6) 

+T-g(a,b). (^(b,&)) (5) 

Using this chain rule, we see that 

x(i,])-G(],f) = *{i,b)-g$,f) + *{i,f) (6) 

where 

m) = L(i, i) • g _1 (i, m) 
*(i,m) = -L(i,i).^(i,m) (7) 

Thus the two fundamental functions of this formalism $, \1/ are closely connected as they 
arise from applying the same operator to the two factors of Q. Defining Yo(i,j) = 
(-t(i,j) + lJ(i,k) 6{i,j))l,*adY 1 {i,j) = t(i,j) (7«+70'))-H(i,j) J(z,k) (7«+7(k)), 
also denote the Fermi gas Greens function 

So% f) = {-(d n - A*)l - Vi)8{i, f) - Y (i, /)}. (8) 
Collecting everything, the exact EOM can now be written neatly as 

{go % j) - A Y x {i J) - A j)} ■ g(j, f) • /i(f, /) 

= tf(i,/)(l-A 7 (<))+A*(<,/). (9) 

We have introduced the parameter A above, with < A < 1, in order to provide an adiabatic 
path between the Fermi gas at A = and the ECFL at A = 1, and also an iterative scheme 
in powers of A connecting the two endpoints. 

We now choose the hitherto undetermined function \x as: 

/i(z,/) = 5(i,/)(l-A7«) + A*(z,/), (10) 
so that Eq. ^ reduces to a canonical FL type equation: 

{go^J) - A n(<,3) - A Hi,])} ■ g(J,f) = 5(i,f). (11) 

Notice that the right hand side has a pure 5 function as in a canonical Fermi liquid type 
theory. To summarize, the EOM Eq. ([T]) under the decomposition Eq. ^ leads to Eq. 
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In turn this splits exactly into two coupled sets of equations Eq. ([7p, Eq. (10) and Eq. (11) 
for the two factors g and /i. Note that the entire procedure is exact, we write explicit forms 
of these equations below and then introduce approximate methods to solve them . 
Inverting we find Dyson's equation for the auxiliary FL Greens function: 



-l, 



i,m) 



{g \i,m) - A Yi0',m) - A $(z,m)}. (12) 



Taking functional derivatives of Eq. (10) and Eq. (12) w.r.t. V, and comparing with Eq. Q 
and Eq. <^ we generate two parallel hierarchies of equations for g and /i that form the core of 
this formalism. The hierarchy for g is essentially autonomous and drives that for /i. Starting 
with the Fermi gas at 0(A°), an iterative process similar to the skeleton graph expansion 
of L-W[l] can be built up, such that terms of 0(A n ) arise from differentiating lower order 
terms of 0(A n_1 ). Systematic approximations may thus be arranged to include all terms of 
0(X n ) for various n[8\. The number of particles is given by \n{i) = g(i, i) • i), and with 

ZC^(a,6;c) = ^ ( S 6) , (13) 
the equations to solve simultaneously are Eq. ([7]), Eq. (12) and Eq. (10). The density 



and spin density response functions (I-Fl,I-F-7) can be found from differentiating Q i.e. 

T^(p, 9 ;r) = ^4 {Gw(p,q)}- 

Zero source limit in Fourier space: When we turn off the source V, the various 
matrix function Q,g,fi become spin diagonal and translation invariant so we can Fourier 
transform these conveniently. We note the basic result expressing Q as a simple product of 
two functions in k space: 
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G{k)= s (k)iM(k), »(k) = i-\- + \*(k) 

(k) = icok + A* - e k {\ - A n) - A (14) 



where e k is the Fourier transform of the hopping matrix —t(i,j), and an uninteresting 
constant term is absorbed in \x here and below. 

Here g plays the role of an underlying auxiliary FL with a self energy <3>, and ^ acts as an 
extra spectral weight that vanishes at high frequency, leaving the exact weight 1 — | valid 
for a projected electron (as in I) for A = 1. Denoting J2k ~^ W~p ^2iui k k w hh N s sites, the 
particle number sum rule is J2 k fi(k)g(k) = |, i.e. 

£ = £g(A0 + A £(*(*) -£)g(*)- (15) 

k k 



In this formalism, at k ~ kp, x = that is relevant to the L-W sum rule, the 3?e g(k) 
dominates 5ie Q(k, 0) (since 9fte \l/(/c,0) is smooth through the FS). Requiring consistency 
with the L-W theorem forces us to pin any sign change of Jie g(k, 0) to the free case, whereby 
we impose a second level sum rule 



£>(*)g(*) = T' and = ?■ (16) 



This can be viewed as a splitting of the usual number sum rule Eq. ( 15 ) [9]. With E(pi,p 2 ) = 
(e Pl + e P2 + | J(0) + | J(pi - p 2 )J we find 

= ^ E(k,p)g(p) A^(p,k) 
p 

= ]T g(p)U^{p,k) (17) 

p 

and the spin labels are from I with the usual significance A* ) = A^ 2 ) — A^ 3 ^ = — |AW. 

Next we introduce the spectral representation of various functions Q that vanish at in- 
finity: Q(iug) = j°° dx and Pq(x) = — - ^sm Q(x + i0 + ), with x + = x + z0 + . The 
Matsubara frequency uq is Fermionic (Bosonic) if Q is Fermionic (Bosonic). Proceeding 
further, at any order in A, the two hierarchies give us coupled equations for the spectral 
densities of the physical particles pg(k, x) as well as the underlying Fermi liquid p s (k, x), in 
terms of the two objects p$(k,x) and p^(k,x) and their Hilbert transforms. The Lehmann 
representation implies that pg(k,x) is positive at all k,x. In making approximations, this 
important and challenging constraint must be kept in mind. 

Solution of g _1 and p to order 0(A) 2 : We next discuss a systematic expansion in 



powers of A [8] , obtained by taking functional derivatives of Eq. ( 10 ) and Eq. ( 12 ) to generate 



expressions for the vertices given the Greens functions via A ~ — ^jg 1 and U ~ -j^p. To 
lowest order in A , the bare vertex A^ a ^ = — 1, this term is absorbed in a renormalization of 



the band dispersion to in Eq. (14) [10J , and the remaining term denoted by $(/c). To 



this order = 0. Proceeding to the next non trivial order in A, by taking the functional 



derivative of Eq. (10) and Eq. (IT2J) we find after a brief calculation: 



V(k) = -2A J2 E i k ,p) Sip) S(q) g(q + P~k) 

p,q 

$(fc) = -2A E(k,p) {E(p,k) + E{q + p - k,p)) 

P,Q 

g(p) g(g) s(q+p- k). (is) 



From Eq. (14) we note that these expressions Eq. (18) lead to a calculation of g -1 and 
p correct upto 0(A 2 ). Frequency dependent corrections arise only to second order in A, 
which is analogous to the structure of the canonical many body theory within the skeleton 
graph expansion. We may now set A = 1 and study the resulting theory as the first step in 
exploring this formalism. 

Denote f(x) = ( exp/ 3 x ) +1 as the Fermi distribution functions and f(x) = 1 — f(x) , and 
denote the usual Fermi factors from second order theory 

W = {f(u)f(w)f(v) + f(v)f(u)f(w)}5(u + w-v-x), 

a function of the frequencies u, v, w, x, and 



y= / W p g (q,w)p s (p,u)p s (q + p- k,v), (19) 

J u,v,w 

a function of k,p,q and x. We may then write the spectral functions corresponding to 



Eq. (18) 



p$$,x) = 2j2E(k,Pl(E(p,k) + E(q + p-k,p))y 

p,q 

P *(k,x) = 2 E(k,p)y. 



(20) 



p,q 



The functions appearing in Eq. (20) are familiar from Fermi liquids [U [5], and encode the 
usual phase space constraints of that theory. This leads to the low temperatures behaviour 
~ max{x 2 , (7t/cbT) 2 }, for both objects 3m ty(k,x,T) and ^sm $>(k,x,T). The real parts of 
these objects are smooth through the Fermi surface, as one expects from the real part of the 



self energy in a FL, and hence motivates the second level sum rule Eq. (16) 



From Eq. (14) we write the exact expression for the physical spectral function pg\ 



pg(k,x) = p g (k : 



1 - 



n 



r)(k, x) 



1 } A(k,x 

where ^ = £k — and the important energy scale A(k, x) and the term rj is defined as: 



(21) 



A(£,x) = 

p-^(k, x) 
rj(k,x) = $t,e^(k,x~ { 



A(k,x 



-Me$(k,x^ 



(22) 
(23) 
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The sign of the energy scale A in Eq. (22) is expected to be positive from Eq. (20). The 



dimensionless term 7] augments the spectral weight at the Fermi level. The equations nec- 
essary to solve the theory to 0(A 2 ) may be summarized as Eq. (14), Eq. (16), Eq. (18) and 



Ref. (JTUl) giving rise to the spectral function Eq. (21). These require further numerical work 
that is underway, it leads to spectral functions in 2 and 3 dimensions that will be published 
separately. However it also provides a very interesting insight about the theory in high 
dimensions that is pursued analytically next. 

Solution in high dimensions: In sufficiently high dimensions, we show next that the 
dimensionless term rj vanishes identically leading to a great simplification. For sufficiently 



high dimensions we can ignore the momentum dependence of y in Eq. (19) and assume 
p$>(k,x) ~ C$ a(x), and p$(k,x) ~ Cq, <r(x), as functions of frequency only. Here a(x) 
extends over energy range u c ~ 0(2W), and C$ has dimensions of inverse energy and is 
positive due to p$ . Its Hilbert transform is called h(x) = V J dy We use an analytically 
tractable Fermi liquid model[UJ with r = irksT, where we set: 



a(x) = {x 2 + r 2 }e 



-C*{X 2 +T 2 }/U C 



(24) 



The peak value of C$ a(x) is of 0(1) and independent of C$ [12J. The other constant is 
dimensionless and negative. To complete the model, we note that the real parts are given 
in terms of h(x) as 9fte $(x + ) = Cq> h(x) and 3fte\l/(x + ) = h(x). With this choice the 
auxiliary spectral weight rj(k,x) vanishes identically in Eq. (23). With T(x) = 7rC$ ct{x) 
and e(£,x) 



[x 



£ — C$ h(x)) we may write p g (£, 



x 



i 



r(x) 



and 3?e g(^, x) 



Denoting (Q(O)c = / d£ Nb(£)Q(£), where Nb(£) is the band density of states 



per spin, the chemical potential is fixed using | 



IZc dx f( x )(P8& x ))t- 



The energy parameter A(k,x) in Eq. (22) is a constant. We scale out a factor to define 

,2 



n 



A = —A{k,x) 



T C 9 ' 



(25) 



The physically observable electronic spectral function reads 

r ( ,) ({!-« + (*){£ 
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(26) 



T 2 (x) + e 2 (^x) 

Here the condition (/)+ = max(0, /), is inserted in the ECFL factor to guarantee the 
positivity of the spectral function for x 2> £p2]. We can determine Aq directly from the 
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second level sum rule Eq. (16): 



/oo 
dxf(x) MZ,x){Z-x})t. (27) 
-oo 

Thus 2/n x Ao is the average inelasticity ||(£ — x)\ \ of the FL Greens function over the entire 
occupied band. It vanishes if p g were a pure delta function, as in a Fermi gas, but is non zero 



in a Fermi liquid. The linear energy term in Eq. (26) thus fundamentally arises to provide 
the extra density to pg, compensating the spectral depletion due to the first factor 1 — | ( 
originating in the non canonical nature of the projected electrons (I)). 

In the numerical solution of the model, we can vary the shapes of the spectra from 
sharp to broad by controlling the energy scale Ao via the parameters C$ and uo in the FL 
function a(x). For illustration we neglect the distinction between the band energy and the 
renormalized choose a flat band density of states per spin po(e) = ^Q(W 2 — e 2 ) hence 
the band width is 2W. Choose C$ = 1 W = 10 4 K pi], this gives A ~ 600K in the cases 



studied. The spectral shapes from Eq. (26) have a characteristic left skew that is visible in 
Fig. ([I]), and also in many experimental spectra in high T c systems. The marginal Fermi 
liquid hypothesis [15] assumes a linear correction to the spectral function, but is symmetric 



about the Fermi energy, i.e. of the form |£ — x\ instead of the term in Eq. (26). 



From Eq. (27) a fascinating duality emerges between the FL and the ECFL[J6]. When 
the FL is overall sharp such that Aq is small, the ECFL is significantly broadened. This 



happens since in the ECFL factor in Eq. (26), the coefficient of £ — x becomes large and 



dominates the 1 — § contribution. The function A(k) in Eq. (22) could vanish at points in 
k space in the full theory (without the assumption of k independence). At those points the 
ECFL spectra would lose all coherence by this duality. A loss of coherence would inevitably 
suggest a (false) pseudo gap, if our current viewpoint were unavailable. The linear term also 
leads to a sloping term in the local density of states of the ECFL that the STM technique 
would probe, although its magnitude and sign are less reliably computed- depending as they 
do on the high energy scales W and ojq. In conclusion, we have presented essential ideas 
underlying the theory of extremely correlated Fermi liquids. We have shown that an explicit 
low order solution is very promising in the context of explaining the photoemission spectra 
of the cuprate materials. 

Detailed numerics and comparison with experiments are currently underway. This work 
was supported by DOE under Grant No. FG02-06ER46319. 



9 




-0.30-0.25 -0.20-0.15 -0.10-0.05 0.00 0.05 



FIG. 1: The density n = .85 and ojq = 0.25. From left to right pg(x) for energies (in units of 
W) £ = -0.3, -0.2, -0.1,0., 0.05 for both the FL (dashed) and the ECFL(solid) theories. Inset 
in (A): provides an enlarged view of the £ = —0.1 plots after inversion, and displays the left-skew 
asymmetry of the ECFL spectrum relative to the FL. Inset in (B) shows the DC resistivity p(T) 
within a bubble approximation as a function of T for the FL (blue) and the ECFL (red). Due to 
spectral redistribution, the ECFL reaches linear T behaviour at a lower T than the FL. 
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